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ABSTRACT 

Context. Variations in the dust emissivity are critical for gas mass determinations derived from far-infrared observations, but also for 
separating dust foreground emission from the Cosmic Microwave Background (CMB). Hi-GAL observations allow us for the first 
time to study the dust emissivity variations in the inner regions of the Galactic plane at resolution below 1°. 

Aims. We present maps of the emissivity spectral index derived from the combined Herschel PACS 160 yum, SPIRE 250 //m, 350 
fjm, and 500 fjm data, and the IRIS 100 /jm data, and we analyze the spatial variations of the spectral index as a function of dust 
temperature and wavelength in the two Science Demonstration Phase Hi-GAL fields, centered at 1=30° and 1=59°. 
Methods. Applying two different methods, we determine both dust temperature and emissivity spectral index between 100 and 500 
//m, at an angular resolution (9) of 4'. 

Results. Combining both fields, the results show variations of the emissivity spectral index in the range 1.8-2.6 for temperatures 
between 14 and 23 K. The median values of the spectral index are similar in both fields, i.e. 2.3 in the range 100-500 yum, while the 
median dust temperatures are equal to 19.1 K and 16.0 K in the 1=30° and 1=59° field, respectively. Statistically, we do not see any 
significant deviations in the spectra from a power law emissivity between 100 and 500 /urn. We confirm the existence of an inverse 
correlation between the emissivity spectral index and dust temperature, found in previous analyses. 
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OS ■ 1. Introduction dPesert et al.L I2008I) found evidence of an inverse relationship 
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I V (A) oc eo — B V (A, T4), (1) 1 in the temperature range 7-20 K, with a behavior similar to 

' that derived from ARCHEOPS data. These variations of B with 

where l y (A) is the specific intensity or brightness, B v is the dust temperature could be owing to intrinsic properties of amor- 

Planck function, T d is the dust temperature, e is the emissivity phous dust grains, as proposed by |Meny et alj (|2007[), but the 

at wavelength A (h and B is the emissivity spectral index. " impact of the noise and the temperature mixing along the line of 

Previous observations in the submillimeter (submm), at ar- sight (LOS) must be carefully taken into consideration, before 

cminute resolution, su ch as those by the balloon-borne ex- a physical interpretation can be made. Variations of the apparent 

periments PRONAOS dDupac et all [2003) and ARCHEOPS spectral index i n a sample of co res in the Galactic plane have 

been observed dHill et all |2006|) using SCUBA and SIMBA. 

* Herschel is an ESA space observatory with science instruments However, because of the restricted wavelength range of their ob- 

provided by European-led Principal investigator consortia and with im- servations, these authors could not attribute them to temperature 
portant participation from NASA. 
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Fig. 1. Left: spectral index in the 100 to 500 fim range, versus dust temperature for each SDP field, combining IRIS 100 fim with 
PACS 160 and SPIRE 250, 350 and 500 fim data, using the least-square fit method in red and the MCMC method in blue. Middle 
and right: spatial distribution of the emissivity spectral index and reduced^- 2 , respectively, derived from the least-square fit method, 
for each SDP field. 



variations. Moreover, none of the previous analyses allowed the 
study of the Tj - /3 correlation at arcminute scales in the inner 
regions of the Galactic plane. Thanks to the Herschel Hi-GAL 
data (9 = 37", 2°x 2° maps, centered in th e Galactic plane at 
1=30° and 1=59°, see lMolinari et all 1201 Oal b), we can now ex- 
tend for the first time this type of analysis to the inner Galactic 
plane, and on continuous fields of diffus e emission. An over view 
of the Herschel instruments is given in Pil bratt et al.1 (l2010h . In 
this paper we analyze emissivity variations between 100 and 500 
fim. We then derive both the dust temperature and spectral in- 
dex from Eq. (Q]i, assuming a single temperature along the LOS. 
This assumption is relatively correct for the 1=59° field, where 
there is less contamination along the LOS in comparison with the 
1=30° field. Indeed, the LOS toward the inner 1=30° field crosses 
several regions with a presumably wider range of temperatures. 
An accurate temperature determination requires sampling both 
sides of the emission peak. The 70 fim brightness is generally 
contaminated by out-of-equilibrium emission from very small 
grains (VSGs) and cannot be used to derive thermal dust tem- 
perature without an accurate subtractio n of this component. Thi s 
contribution has been investigated by Compiegne et al] (12.010b . 
Therefore we include the IRAS data at 100 fim in our analysis, 
for which we estimate an average VSG contamination of less 
than 10%. This implies degrading the original resolution of the 
Hi-GAL data by adopting the IRAS angular resolution of 4'. We 
note, as a consequence, an averaging effect within the resolu- 
tion element, which also needs to be taken into account in the 
interpretation of the results of the present analysis. 

One of the key aspects in the analysis of CMB is the ability 
to separate its emission from the other astrophysical foregrounds 
(including thermal dust) through multifrequency observations. 
While it is not easy to achieve successful component separation 
over the Galactic plane, determining the dust spectral index vari- 
ations across the sky can effi ciently help in reducin g the number 
of unknowns in the problem (Ricciardi et al .11201 Oh . 

2. Data 

We use the RO MAGAL Hersc h el PA CS and SPIRE 
maps described in iTraficante et al.l (1201 Oh . combined with 



the IRIS (Improved Reprocessin g of t he IRAS Survey, 
see Miville-Deschenes & LagacheL [2005) 100 fim data. 
Multiplicative factors (0.78,1.02, 1.05 and 0.94 at 160, 250, 
350 and 500 ^ m , resp e ctively) have b een applied to the data 
dPoglitsch et all l2010t iGriffin et all l2O10t ISwinvard et ail 
2010). An absolute calibration accuracy of 20% and 15% for 
PACS a nd SPIRE has been adopted. We applied the offsets 
given in iBernard et al.1 {2010) All maps have been convolved 
with a Gaussian kernel, with a FWHM of 4', i.e. equal to the 
IRIS 100 fim angular resolution. In addition, the maps have 
been rebinned on a uniform grid with a pixel size of 1.65'. The 
IRIS 100 nm calibration uncertainty is taken to be 13.5% (see 
Miville-Deschenes & Lagache, 2005). 



3. Variations of the emissivity spectral index with 
temperature 

We use two different methods for the simultaneous derivation of 
the dust temperature and the emissivity spectral index. Both al- 
low us to fit the data with a modified blackbody function (see 
Eq. [TJ. In the first method we perform an y 2 minimizatio n, ap- 
plying the same code as described in [Dupacet al. (2001), us- 
ing the IDL least-square fit curve fit function to deduce the Td 
and fi parameters, as well as their lcr uncertainties. The second 
method estimates the best-fit parameters by looking for the max- 
imum like lihood, using a Monte C arlo Markov Chain (MCMC) 
algorithm (Lewis & Bridle, 2002) and represents a bayesian ap- 
proach to parameter estimation. The posterior distribution for 
the parameters is sampled by using the Metropolis-Hastings al- 
gorithm, and a maximum likelihood estimate is derived jointly 
for /3 and Td- We chose a wide flat a priori probablity density 
of the parameters, i.e. K< Tj <60 K and -1< (5 <5, in order 
not to constrain the fit results. A benefit of this method is that it 
also recovers the joint posterior distribution of the two estimated 
parameters. We make use of this feature to explore the T^-f3 cor- 
relation as explained below. 

We fit the spectral energy distribution (SED) between 100 
and 500 fim for each pixel of the maps, using both methods. The 



see their Table 1 
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Fig. 2. Representative SEDs combining IRIS 100 /jm with PACS 160 /jm and SPIRE 250, 350 and 500 yum, modeled by the gray- 
body law with Td and ft derived from the least-square fit method between 100 and 500 fim (left and central panels), and between 
100 and 350 //m (right panel). 



color correction is computed iteratively in all channels. Results 
of the Td-/3 determination for each field are shown in Fig. Q] 
Only pixels with a surface brightness at 500 //m higher than 
105 MJy/sr and 50 MJy/sr for the 1=30° and 1=59° fields, re- 
spectively, have been considered. Both methods show a clear 
anti-correlation Tr/3 in each field. These two independent meth- 
ods also highlight a good agreement. The ft values range from 
^ 2.6 to ^ 1.9 for temperatures between 17.5-23 K for the 
1=30° field, and from 2.7 to =* 1.8, between 14 and 21 K, 
for the 1=59° field. The assumption made of a single tempera- 
ture along the LOS may hold in the 1=59° field, but is certainly 
too simplistic for the 1=30° field, since tempera ture variations 
along the LOS are larger in the inner field (see iBernard et all 
l2010t) . The analysis of the 1=30° field would certainly require 
an accurate combination of temperatures and column densities. 
However, our methods do not allow the simultaneous fitting of a 
complex mixture of spectra. Even if the 1=30° field is potentially 
characterized by several grain temperatures along the LOS (as- 
sociat ed to Sagit t arius, Scuttum-Crux and Perseus, as described 
in IBernard et all 1201 Oh . and therefore cou l d induce a spuriou s 
T d -ft inverse correlation dMasi et all 119951: IShettv et all l2009h . 
the 1=59° field, less contaminated by various dust mixing ef- 
fects, still highlights a pronounced behavior in the T^-ft param- 
eter space. Both fields, with a different mixture of temperatures, 
present a T^-ft anti-correlation. We therefore interpret this as an 
indication that the mixture of temperatures is probably not the 
dominant responsible effect for the observed anti-correlation. 

The spatial distribution of ft is shown in Fig. Q] Warmer 
regions, which most often correspond to bright regions, show 
lower beta values. However, the 1=59° field presents a gradient 
along the East- West direction, which is not observed in the PACS 
or SPIRE surface brightness maps. The stripes visible in the 
maps are associated to low-level residual stripes in the IRIS data. 
In the 1=59° field, ft values higher than 2.6 located at 1=57°48', 
b=-0°20' correspond to one of the less active star-forming re- 
gions, although the PACS 160 data show an artifact at this loca- 
tion, disclosed by a high x 2 (see Fig. Q]). In the same field, the 
lowest ft values (< 1.8) are found for an HII region, with a bright 
source nearby at 1=59°38', b=+0°38'. Representative examples 
of SEDs with a gray -body fit are presented in Fig. [2] illustrating 
the agreement between the data and the model. 

The two parameters Tj-/3 are degenera te in the parameter 
space, creating a spurious inverse relation ( Sh ettv et all [2009) 
that has to be considered while investigating if an intrinsic 
physical correlation indeed exists. In order to properly estimate 
the relationship between Td and ft, we have to take into account 
the correlation that is shown by the spectral shape in Eq. ([TJ 
and caused by calibration uncertainties. One way to include the 



effect of the degeneracy is to estimate this effect not only for the 
Td-/3 best-fit values but for points included in the 68% contours 
of the two-dimensional posterior probability obtained with the 
MCMC method. Within the 68% contours the points are not 
uniformly distributed, but their density is higher close to the 
maximum of the probability, increasing their weight, and lower 
close to the boundaries. Assuming the model ft = A(Td/2QK)~ a , 
between 100 and 500 //m in both fields, the fit is then performed 
on all the T^-ft couples inside a posteriori probability contour of 
each pixel, chosen randomly, one for each pixel and weighted as 
described. We then take a point within the contour for each pixel 
of the map and fit the trend to estimate the A and a parameters. 
We repeat this procedure until we obtain a good sampling of 
the shape of the distribution of posterior probabilities. This 
method has been succ essfully tested on BOOMERanG data 
dVeneziani et all 1201 Oh and allows us to include the effects of 
degeneracy in the error bars. The resulting A and a ditributions 
are well approximated by Gaussian functions, even if a slight 
asymmetry is present. The center of the Gaussian is the more 
likely value, while the error is estimated by marginalizing over 
the other parameter. The best-fit values so obtained are 



,j >-1.32±0.04 

ft = (2.19 + 0.25) x ,/ = 30° (2) 



,j >-1.33±0.04 

ft = (1.71+ 0.22) x U ,/ = 59°, (3) 



Figure [3] shows the results of the fitting method over the prob- 
ability 68% contours with only a few pixels out of ~ 3500 
contours for clarity. The best-fit Hi-GAL Tr/3 relation is plot- 
ted, as well as the ARCHEOPS, BOOMERanG, and PRONAOS 
ones for comparison. Our results highlight a steeper Tr/3 law 
in the inner Galactic plane in both fields than in external 
regions and high latitude clouds observed with ARCHEOPS 
and BOOMEranG. Whereas the wavelength range covered by 
PRONAOS is almost the same as that of Herschel, the T&-ft re- 
lation (ft - 1/(6 + cjTd)) is flatter, with ft values going from 2.0 
to 1.7, in the range 13-26 K. These results could indicate dif- 
ferent dust properties in the inner Galactic plane, as compared 
to the outer galaxy observed with ARCHEOPS, the solar neigh- 
borhood observed with PRONAOS, and high latitude clouds ob- 
served with BOOMEranG. We therefore think that these differ- 
ent Td-/3 relations could be caused by changes in dust properties 
with the environments. 
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Fig. 3. Td-/3 two-dimensional 68% contour posterior probabili- 
ties derived from the MCMC method, shown for ^15-25 pixels 
for the two SDP fields. The overplotted lines correspond to Hi- 
GAL (solid red and blue lines for the 1=30° and 1=59° fields, 
respectively), ARCHEOPS (dashed line), BOOMERanG (line 
with circles) and PRONAOS (line with squares) best-fits, respec- 
tively. The Td-y6 relationship is estimated over all points within 
the contours. The gray points correspond to the T^-ft data points 
derived from the MCMC method (see Fig. 1, left panels). 



4. Variations of the emissivity spectral index with 
wavelength 

In the previous section we have found evide nce of spectral inde x 
variations with temperature. As predicted by Meny et al. (2007), 
a flattening of the emissivity spectrum (and so a lower ft) is ex- 
pected for wavelengths longer than 500 fim as a consequence 
of the internal structure of the grains. Recently, iGordon et alj 
(2010) highlighted an emission excess at 500 fim in the Large 
Magellanic Cloud that could be of the same origin. However, this 
excess is of about 10%, which is lower than the SPIRE calibra- 
tion uncertainties of 15%. The Herschel data do not allow us to 
study emissivity variations at longer wavelengths than 500 fim. 

Although the SEDs in each field (left and central panel of 
Fig- EJ d° not provide evidence for any departures from a power 
law, we analyze the possibility of an emissivity excess at 500 
fim. To properly perform this study, the temperature estimates 
need to be done independently of the 500 fim data. Therefore 
we have redone the Tj and ft determination with the least-square 
fit method using only data between 100 and 350 fim, and we 
compare the emission at 500 fim with the value predicted by ex- 
trapolating the emissivity power law between 100 and 350 fim. 
The predicted values are systematically lower than the data at 
500 fim in both fields, by 16% in 1=30° and 13% in 1=59°, 
which could favor the hypothesis of an emissivity excess at 
this wavelength. These values are not significant compared to 
the uncertainties on the SPIRE data. Moreover, we suspect that 
the PACS 160 um data suffer from calibration issues, as dis- 
cussed in iBernard et al.l d2010h . The calibration uncertainty at 
160 fim does not impact the results of the global analysis de- 



scribed above, but here the limitation of the wavelength range 
(100-350 fim) removes some constraints on the fit, resulting in 
an overestimate with respect to the data at 250 fim, and an un- 
derestimate at 500 fim (see Fig. [2] right panel). However, if we 
slightly shift down the 160 fim data points, the apparent under- 
estimate at 500 fim would disappear. For this reason, we do 
not claim that this excess is real and we think instead that the 
emissivity spectra are relatively constant between 100 and 500 
fim with similar median ft values in both fields, i.e. 2.3, even if 
the median temperature is statistically lower in the 1=59° field 
(T d = 19.1 K for 1=30° and T d 16.0 K for 1=59°). 

5. Conclusions 

We investigated variations in the spectral index of the dust emis- 
sivity, with temperature and wavelength, in the inner Galactic 
plane, using the new Herschel observations in two Hi-GAL 
fields, centered at 1=30° and 1=59°, acquired during the Herschel 
Science Demonstration phase, combined with the IRIS 100 fim 
data. We fitted the SEDs for each pixel of the two fields with two 
independent methods, deriving simultaneously the emissivity 
spectral index and the dust temperature by adjusting a modified 
blackbody function to the data. The results are similar with both 
methods. Using a Monte Carlo Markov Chain algorithm method 
we computed the 68% likelihood contours for each point. We 
find a Td-ft inverse correlation, with the local variation going 
from 1.8 to 2.6 for temperatures between 14 and 23 K, shown 
for the first time in the inner Galactic plane. Moreover, our re- 
sults indicate a different trend with respect to previous investi- 
gations based on BOOMERanG, ARCHEOPS, and PRONOAS 
data, probably because of different dust properties in the inner 
Galactic plane compared to other environments. An extensive 
follow-up analysis will be required to take into consideration 
the mixture of temperatures along the line of sight. The median 
value of ft is similar in both fields, equal to 2.3, slightly higher 
than the usual reference value of 2. We do not favor the hypoth- 
esis of an emissivity excess at 500 fim, as suggested for the 
Large Magellanic Cloud. A complementary forthcoming study 
will combine Planck with Herschel data to investigate possible 
changes in the dust emissivity spectrum at wavelengths larger 
than 500 fim, as reported in previous studies of dust emission in 
our Galaxy. 
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